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1. Introduction 


New physics beyond the Standard Model (BSM) in the framework of some new strongly- 
interacting gauge theory with a composite Higgs mechanism is an attractive BSM scenario with 
related lattice work reviewed recently [1]. Of course we hear voices that pursuing the composite 
Higgs scenario is overtaken by recent findings at the LHC. After all, a light Higgs-like scalar was 
found, consistent with SM predictions, and composite states have not been found below the TeV 
scale. In contrast, the voices continue, strongly coupled BSM gauge theories are Higgs-less with 
resonances predicted below the TeV scale, just like in the original technicolor idea [2, 3]. The 
facts do not seem to support this skeptical view which originates from naively scaled properties 
of Quantum Chromodynamics (QCD) to the TeV region. Related old technicolor guessing games 
were lacking any credible predictive power close to the conformal window where gauge theories 
are nearly scale invariant, in sharp contrast to QCD which is not. In fact, there is no evidence 
that compositeness and a light Higgs scalar are incompatible. Recent developments are hinting 
compatibility, like in near-conformal gauge theories where a light composite scalar could perhaps 
emerge on the Electroweak scale with a resonance spectrum far separated above the TeV scale, 
perhaps within the reach of Run 2 at the LHC. 

Work in progress is reported here on a particularly interesting gauge theory with a fermion 
doublet in the two-index symmetric (sextet) representation of the SU(3) color gauge group [4, 5, 
6, 7]. Our observations suggest that the model is very close to the conformal window and a light 
composite scalar appears to emerge with 0++ quantum numbers. From chiral symmetry breaking 
we find fhree massless Goldsfone pions in fhe specfrum. Wifh Elecfroweak inferacfions fumed on, 
fhe model exhibifs fhe simples! composife Higgs mechanism and leaves open fhe possibilify of 
a lighf composife scalar sfafe wifh quanfum numbers of fhe Higgs imposfor emerging as perhaps 
fhe pseudo-Goldsfone dilafon-like sfafe from sponfaneous symmefry breaking of scale invariance. 
Even if scale symmefry breaking is enfangled wifh ;^SB wifhouf dilafon inferprefafion, a lighf 
Higgs-like scalar sfafe can emerge from fhe new gauge force close fo fhe conformal window. The 
main goal of our Higgs projecf is fo invesfigafe fhese imporfanf problems wifh ah inifio laffice 
simulafions of fhe sexfef model. 

In Secfion 2 we infroduce fhe Electroweak embedding of the strongly coupled sextet gauge 
theory, comment on the intriguing features of the lowest baryon state in the minimal sextet model 
and its extensions, and describe the new data sets developed since our last report [8, 9]. In section 3 
we investigate Goldstone dynamics and Electroweak scale setting from chiral symmetry breaking 
as premier ingredients of the composite Higgs mechanism. We also analyze cutoff dependent taste 
breaking effects in the non-Goldstone pion-like spectrum of staggered fermions. In Section 4 we 
test the GMOR relation from the spectrum of the Dirac operator and the related chiral condensate. 
In Section 5 we present new results on the light 0++ scalar and outline future plans. In Section 
6 we begin to develop and test mixed action based improved analysis of ;i;SB with new run plans 
at fixed fopology fo cross over from fhe p-regime fo fhe epsilon-regime of ;^SB. In Secfion 7 we 
probe fhe scale-dependenf running coupling from fhe perturbative UV scale fo fhe scale of chiral 
symmefry breaking. 
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2. Electroweak embedding and computational framework 


Quantum numbers and Electroweak symmetry breaking pattern 

The two fermion flavors of the model transform in the eomplex two-index symmetrie (sex¬ 
tet) representation of the SU(3) eolor gauge group whieh implies SU(2 )l x SU(2)r x U(1) flavor 
symmetry for the gauge foree of the theory. The fermions are assembled into a left-handed weak 
isospin doublet qL and two right-handed weak isospin singlets qr with 






u 


(if) 

R ’ 



( 2 . 1 ) 


where i,j = 1,2,3 label the eolor indiees of the symmetrie tensor elements. The ABJ anomaly 
would spoil the renormalizability of the gauge theory so the fermion gauge eoupling must not in- 
troduee anomalous Ward indentities. This requires the traee tr({T®(R),T'’(R)}T‘^(R)) to vanish 
for fermion flavor group representation R with representation matrix T^(R). In the sextet model, 
the fermions are either doublets or singlets under the flavor group SU(2). The matrix T® will 
be either the Pauli matrix or the U(l) hypereharge Y. Sinee the SU(2) group is anomaly-free, 
tr({T‘, = 25‘jtr(T'^) = 0, it is easy to see that tr(Y) = tr(Y^) = 0 are the two anomaly eondi- 

tion to satisfy. The absenee of gauge anomalies thus requires a traeeless weak hypereharge operator 
for the left-handed fermion doublet. Choosing the generator of the U(l) group as Y = 2(Q —T 3 ) 
for the weak hypereharge Y, with eharge Q and the third eomponent T 3 of weak isospin, requires 
half-unit of eharge for the left-handed fermions whieh form a weak isospin doublet. 
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The right-handed u-fermion will have hypereharge Y = 1 and the right-handed d-fermion will have 
hypereharge Y = — 1 with eonsistent eharge assignments 
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satisfying the tr(Y^) = 0 anomaly eondition. The ehiral SU(2 )l x SU(2)r symmetry of the the¬ 
ory is dynamieally broken to the diagonal veetor SU(2)v subgroup. This ;^SB is responsible for 
breaking SU(2)w x U(1 )y to U(l)em- The residual SU(2)v symmetry is the eleetroweak analog of 
isospin and the approximate weak isospin invarianee of the eleetroweak foree and the new gauge 
foree ensure that the p-parameter is approximately one. There is also an exaet U(l) symmetry in 
the theory whieh is unaffeeted by ehiral symmetry breaking proteeting the eonserved baryon num¬ 
ber. Baryon number and eharge eonservation are keeping the lightest baryon stable against the new 
gauge foree and weak deeays. 


Baryon construction in the sextet model 

The eharge assignment has intriguing implieations for the baryon speetrum of the sextet model. 
The sextet representation of fermions with SU(3) eolor gauge group will impose a symmetrie eolor 
wave funetion for baryon states as three-fermion systems. This is in sharpe eontrast to QCD where 
the eolor wave funetion of the nueleon is antisymmetrie. The non-relativistie limit of the flavor- 
spin-spatial part of the baryon wave funetion will look like the wave funetion of triton in terms of 
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symmetries [10]. It follows from the construction of Eqs. (2.1-2.3) that the lightest baryons form a 
stable isospin doublet of (uud) and (udd) states which carry half-integer charges with opposite sign. 
Detailed properties of the sextet baryons with lowest mass in the 3 TeV range were reported at the 
conference [11]. Dark matter related relic sextet baryon issues are expected to come into focus only 
if the model will deliver a viable composite Higgs mechanism which remains the primary focus of 
our investigations. 

Dark matter 

The lowest stable baryon state in the sextet model is in the fractionally charged massive par¬ 
ticle (FCHAMP) category of dark matter speculations regarding the evolution history of the Uni¬ 
verse [12, 13]. It is difficult and ultimately imperative to estimate the relic abundance of the frac¬ 
tionally charged and stable sextet baryon. Only qualitative arguments can be given that the very 
small relic abundance is likely to escape existing experimental limits and theoretical requiem [13]. 
Sextet baryons and antibaryons are produced in the Electroweak phase transition which is expected 
to be of second order with two fermion flavors in the chiral limit. Charge symmetric baryon and 
antibaryon densities in thermal equilibrium will continue to decrease well below the critical temper¬ 
ature Tc until at freeze-out temperature T* the expansion rate will set the density to its relic abun¬ 
dance level from the solution of the Boltzman equation [14]. The freeze-out temperature and the 
related relic abundance level are very sensitive to the annihilation rate of baryons and antibaryons 
in the sextet model. As a general trend, the stronger the interaction, the longer the particles remain 
in equilibrium (larger x* = Mb/T*) and the fewer survive (~ e^’'* ). Earlier technicolor estimates 
which were based on scaled up QCD calculations of annihilation cross sections are irrelevant for 
the new theory which is close to the conformal window. The expected relic abundance is minuscule 
but more quantitative calculations are required with some sense of urgency to settle this interesting 
question. 

Some surviving sextet baryon asymmetry in the early evolution of the Universe could affect 
the above argument based on a symmetric sextet baryon-antibaryon distribution. B-violation as 
required by the first of the Sakharov conditions can be associated with the non-Abelian anomaly 
of the left-handed fermion current in nontrivial background gauge field configurations outside per¬ 
turbation theory [15]. For the second Sakharov condition, the origin of C and CP violation, if any, 
is left undetermined in the new theory and the well-known effects related to sphaleron dynamics 
could wash out any early sextet baryon asymmetry generated before the Electroweak phase transi¬ 
tion [15]. Explaining dark matter would probably require an extension of the minimal sextet model 
by adding a new lepton doublet with quantum number assignments of the doublets and singlets 
following the QCD pattern. It could make the lowest baryon state neutral but still requires efficient 
C and CP violation of ill-understood origin. As an alternate extension, a third fermion flavor could 
be added which is massive and remains an Electroweak singlet. 

Algorithms, codes, and run parameter sets 

We use the tree-level Symanzik-improved gauge action for all simulations reported here. The 
conventional jS = 6 /lattice gauge coupling is defined as the overall factor in front of the well- 
known terms of the Symanzik lattice action. The link variables in the staggered fermion matrix 
are exponentially smeared with two stout steps [16]; the precise definition of the staggered stout 
action was given earlier in [17] and the RHMC algorithm has been deployed in all runs. The 
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fermion flavor doublet requires rooting in the algorithm. For molecular dynamics time evolution 
we apply multiple time scales [18] and the Omelyan integrator [19]. We have highly efficient codes 
running on BG/Q, gpu, and commodity cluster platforms. Our error analysis of hadron masses 
is based on correlated fitting with double jackknife procedure on the covariance matrices [20]. 
The time histories of the fermion condensate, the gauge field energy on the gradient flow, the 
topological charge, and correlators are all used to monitor autocorrelation times in the simulations. 
We have new simulation results at j3 = 3.2 and 3.25 for fermion masses m = 0.002,0.003,0.004 
on 32^ X 64, 40^ x 80, and 48^ x 96 lattice volumes. We also have new runs at jS = 3.15 with 
m = 0.003,0.004,0.006,0.008 and at j3 = 3.30 with m = 0.005,0.006,0.008,0.010 on 32^ x 64 
lattice volumes. 


3. Goldstone spectrum and Electroweak scale setting 


If the chiral SU(2 )l x SU(2)r symmetry of the model is dynamically broken to the diago¬ 
nal vector SU(2)v three associated Goldstone pions facilitate the minimal realization of the Higgs 
mechanism after the Electroweak interactions are turned on. As shown in Figure 1, our results 


Goldstone pion (Rwall pion channel) Goldstone plon {Rwall plon channel) 




F decay constant (Rwall plon channel) 


F decay constant (Rwall pion channel) 




Figure 1; The leading order fits are shown at two values of j8 = 6/gQ for the fermion mass dependence of the Goldstone 
pion from chiral perturbation theory (;tPT) without logarithmic loop corrections. Ignoring taste breaking the fundamen¬ 
tal B parameter of the chiral Lagrangian is given by Cm/ 2 in leading order (LO). The fits of Fj^ include the linear NLO 
analytic term from p^PT. All fits are based on the staggered pion correlator with exact PCAC relation using random wall 
noise vectors. We have similar analysis for and F;t at jS = 3.15 and 3.30 as well. 

are consistent with chiral symmetry breaking exhibiting consistent Goldstone pion behavior under 
fermion mass deformations. The Electroweak scale in finite lattice spacing units is set from the 
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pion decay constant F;;: in the chiral limit with F = 250 GeV in continuum physics notation. The 
preliminary results of Figure 1 represent work in progress with moving parts which include contin¬ 
uing refinement of the fitting procedures on our large new data set, the unfinished analysis of taste 
breaking in staggered chiral perturbation theory, new runs closer to the p-regime of leading chiral 
logarithms, and the influence of a light scalar state on the analysis of chiral perturbation theory. 
In Section 6 we will briefly sketch new directions with crossover from the p-regime to the epsilon 
regime and Random Matrix Theory (RMT) applying mixed actions in the analysis. 


Taste breaking cutoff effects 

Since the determination of the Goldstone decay constant F in the chiral limit is critically im¬ 
portant for the determination of the light scalar mass and the separated resonance spectrum, we 
will briefly describe taste breaking effects which will influence the final outcome of the analysis. 
To illustrate cutoff dependent taste breaking effects, spectra of selected non-Golds tone pion states 


(3=3.20 non-Goldstone scPion spectrum 



m 


p=3.20 non-Goldstone ijPion spectrum 



m 


p=3.25 non-Goldstone scPIon spectrum 


p=3.25 non-Goldstone ijPion spectrum 




Figure 2: LO fits from the analytic mass dependence of the chiral Lagrangian without logarithmic chiral loop correc¬ 
tions are shown for two non-Goldstone pions at two values of j8 = 6/gQ. We have similar fits at jS = 3.30. The scPion is 
degenerate with the iSPion (not shown) in the same SO(4) multiplet [21]. 


are analyzed in Figure 2 with the definition of the relevant correlators and quantum numbers given 
in [22, 23]. In the fermion mass range of our data set the taste breaking pattern is different from 
QCD where the residual A mass shifts of the non-Goldstome pions are equispaced in the chiral limit 
with approximately degenerate SO(4) taste multiplets and with parallel slopes for finite fermion 
mass deformations of Goldstone and non-Goldstone pion states [21]. For example, as part of the 
equispaced split of degenerate SO(4) multiplets in QCD, the observed approximate split Ay ~ 2Asc 
of two multiplets appears to have collapsed in the sextet model from our fitting procedure. 
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The other distinct difference from QCD is the non-parallel slopes which fan out in Goldstone 
and non-Goldstone mass deformations of the pion spectrum as shown in Figure 2. While the A ad¬ 
ditive mass shifts are LO taste breaking effects in the chiral Lagrangian [21, 24], the taste breaking 
slope corrections 5 can plausibly be identified with NLO analytic terms in the chiral analysis [25]. 
The corrected mass relation is ~ ^lo( 1 + where 5 depends on the taste quantum num¬ 
ber of the pion state. Several relations constrain the 5 taste breaking corrections [25] . For example 
the relation = — 5ij immediately implies that the fitted slope of the Goldstone pion must re¬ 
ceive significant taste breaking and cutoff dependent correction in the linear fit of Figure 1 since 
the slopes of the Goldstone pion and ijPion fan out considerably in Figure 2. We can infer from 
the measured slopes the important relation 5ji = — (cij — CM)/(cij -|-cm) to determine the leading 
cutoff correction to the fundamental B parameter of the chiral Lagrangian from the fit parameter 
Cm as 2B = cm/(1 + Stc)- The correction factor 5yt decreases from -0.207 to -0.087 as j3 is varied 
from 3.20 to 3.30 with decreasing lattice spacing. Work on cutoff corrections to the decay con¬ 
stant F are in progress. The small A mass shifts in the chiral limit and the significant fan-out taste 
breaking structure of the slopes led us to generate a new data set below the fermion mass region 
m = 0.002 — 0.006. This new effort crossing over from the p-regime toward the epsilon regime and 
RMT using mixed action based analysis will be outlined in Section 6. 


4. Chiral condensate and GMOR 


The consistency of the fundamental parameters F and B and the direct determination of the 


spectral density p(k,m) 



time history of Q topological charge 



spectral density p(k,m) 




Figure 3: Tests on the spectral density of the staggered Dirac operator are performed on our largest 48^ x 96 lattice 
volumes at j3 = 3.20 with two different fermion masses. Our new Chebyshev expansion is also shown. 
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non-vanishing fermion condensate £ in the chiral limit is tested by the GMOR relation 2BF^ = £ 
where summation over two fermion flavors is included [26]. Access at j3 = 3.20 to the numerical 
estimate 2BF^ = 0.00497 within the range of a few percent error is provided by the analysis of 
data and fits shown in Figure 1 and Figure 2. The slope correction 5^ is a significant factor in the 
analysis. The Banks-Casher relation [27] relates the condensate £ to the spectral density p(A,m) 
of the Dirac operator, 


p(A,m) = y (5(A — Ayt))/F, with lim lim lim p(A,m) = £/(27r), (4.1) 

where the spectral density is determined as the ensemble average over the Dirac eigenvalue density 
in finite volume V. In Eq. (4.1) the condensate £ = —{\j/Y) for two fermion flavors is determined 
by the eigenvalue density. Figure 3 shows a select subset of recent analysis of the spectral density of 
runs with 48^ x 96 lattice volume at j3 = 3.20 for the two lowest fermion masses m = 0.002,0.003. 
The lower left panel is the topological history of the m = 0.002 run on the gradient flow at flow 
time t = 20. The upper left panel shows the spectral density (with the In factor absorbed) for a 
select subset of gauge configurations with topological charge Q = 0 and the upper right panel with 
Q = — 1. The results on the spectral density p(A = 0,m), at the lowest values of A reached for 
the A —5- 0 limit, are not sensitive to the two values of the topological charge tested. Agreement 
with GMOR is remarkably good from the independent few percent level determination of 2BF^ = 
0.00497 as discussed above. The new analysis removes earlier inconsistencies from the sextet 
GMOR relation [23]. Continued work is necessary for a more complete analysis of the systematic 
effects. Comprehensive finite size scaling analysis, the chiral limit m —0 of the spectral density, 
and the scale-dependent renormalization of the condensate remain important unfinished goals. 

The determination of the spectral density from the low eigenvalues of the Dirac operator has 
a limited range and becomes increasingly difficult for larger volumes. In several applications, 
like the anomalous dimension of the mode number density, it is important to determine p (A, m) 
for a large range of A and in large lattice volumes. Recently we developed and tested a new 
stochastic method with random noise vectors which is capable of calculating the entire spectral 
density function and mode number distribution of the Dirac operator with great efficiency [28]. 
The method is based on a high precision finite resolution Chebyshev approximation to the Dirac 
delta function in Eq. (4.1). At any given value of A in the spectrum we can extrapolate to infinite 
polynomial order in the Chebyshev expansion utilizing well-known asymptotic properties of Tn (A) 
Chebyshev polynomials at fixed A in the n —oo limit. The expansion used here is different from 
the mode number approximation introduced earlier at a fixed value of A [29]. 

The lower right panel of Eigure 3 displays a typical result on our largest sextet lattice vol¬ 
ume, averaged over gauge configurations for the full eigenvalue spectrum. When magnifying the 
low infrared part of the spectral density with the same expansion, the two upper panels show the 
convergence rate of the Chebyshev approximation as a function of the the polynomial order when 
compared with the direct diagonalization of the Dirac operator. Polynomial order in the n = 6000 
range is almost indistinguishable from the data and the extrapolation procedure works well from 
lower orders. There is a variety of interesting applications where this method can be further ex¬ 
plored. 
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5. The 0++ light scalar and the resonance spectrum 

The most important goals of our lattice Higgs project are to establish the emergence of the 
light scalar state with 0++ quantum numbers and the resonance spectrum perhaps far separated 
from the light composite scalar. 

The light scalar state 

The fo meson (in QCD terminology) has quantum numbers and acts as the scalar state in 
the sextet model. Close to the conformal window, the fo meson of the sextet model is not expected 
to be similar to its counterpart in QCD. If it turns out to be light, it could replace the elementary 
Higgs particle and pose as the Higgs impostor. Two types of different operators, the fermionic 
one and the gluonic one (0++ glueball), are expected to mix. Such mixing was not included in the 
pilot study [8] but becomes an important goal of our ongoing effort. 

Triplet and singlet masses from 0'^'^ correlators 
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Figure 4: The left panel shows earlier results on 32^ x 64 lattice volumes at jS = 3.20 [8]. The right panel shows some 
representative new result on a large lattice volume using correlated fitting and Principal Component Analysis (PCA). 

A particular flavor-singlet correlator is needed to capture the 0++ scalar state with vacuum 
quantum numbers. It requires connected and disconnected diagrams of fermion loop propaga¬ 
tors on the gauge configurations. The connected diagram corresponds to the non-singlet correla¬ 
tor Cnon-singiet(t)- The Correlator of the disconnected diagram is D(t) at time separation t. The 
fo correlator Csingiet(t) is defined as Csingiet(t) = Cnon-singiet(t)+D(t). The transfer matrix has 
the spectral decomposition of the Csingiet(t) correlator in terms of the sum of all energy levels 
Ei(0++),i = 0, 1,2,... and their parity partners Ej(0^+),j = 0, 1,2,... but at large time separation 
t the lowest states Eo(0++) and Eo(0 dominate. They correspond to mf^ and m,j^^. The rel¬ 
evant non-singlet staggered correlator can be fitted well with non-oscillating ao contribution and 
oscillating TTsc contribution, with the non-Goldstone pion TTsc discussed in Section 3. 

We estimate the connected and disconnected diagrams with stochastic source vectors of fermion 
propagators. To evaluate the disconnected diagram, we need to calculate closed loops of quark 
propagators. We introduce Z 2 noise sources on the lattice where each source is defined on indi¬ 
vidual time-slice to for color a. The scheme can be viewed as a “dilution” scheme which is fully 
diluted in time and color and even/odd diluted in space. The left panel of Eigure 4 shows the earlier 
preliminary results from the pilot study on 32^ x 64 lattice volumes at j8 = 3.20 [8]. Erom our new 



(5=3.20 (with PCA analysis) 48^x96 m=0.002 
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number of blocks = 16 block size = 4 
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analysis a representative example of the scalar effective mass fit is shown at j8 = 3.20 as the right 
panel of Figure 4 on a large 48^ x 96 lattice volume and probes the light scalar deeper toward the 
chiral limit at fermion mass m = 0.002. 

Further work is needed on the light fo scalar with 0++ quantum numbers. The lower left panel 
of Figure 3 indicates large autocorrelation times in the topological history of the RHMC algorithm. 
The effects of slowly changing topology on the fitted mass values of fo are not sufficiently tested 
although sensitivity to the topological charge remains within the statistical accuracy of the runs. 
Fermion mass deformations of the low-lying fo state and the Goldstone pion are expected to be 
entangled which requires the modification of ;^PT. Precise extrapolation to vanishing fermion 
mass in the chiral limit remains a challenging problem in the presence of the light fo state. 

The emerging spectrum 

It is important to investigate the chiral limit of composite hadron states separated from the 
Goldstones and the light scalar by finite mass gaps. The baryon mass gap in the chiral limit provides 
further evidence for ;i;SB with preliminary results reported at this conference [11]. Resonance 
masses of parity partners provide important additional information with split parity masses in the 
chiral limit. This is particularly important for consistency with ;^SB and for a first estimate of the 
S parameter when probing the model via Electroweak precision tests [30]. 

A remarkable spectrum is emerging which is sketched in Figure 5 for illustration only. Al¬ 
though with more work needed to confirm, the sextet model appears to be close to the conformal 
window and due to ;t;SB exhibits the right Goldstone spectrum for the minimal realization of the 
composite Higgs mechanism with a light scalar separated from the associated resonance spectrum 
in the 2 TeV region. Chiral symmetry breaking and a very small beta function, perhaps slowly walk¬ 
ing as hinted by preliminary results in Section 7, are not sufficient to guarantee a light dilaton-like 
state as the natural explanation for the emergence of the light scalar. Consistent with our observa¬ 
tions, a light Higgs-like scalar is still expected to emerge near the conformal window as a composite 
state with 0++ quantum numbers, not necessarily with dilaton interpretation. This scalar state has 
to be light but is not required to match exactly the observed 126 GeV mass. The light scalar from 
composite strong dynamics gets lighter from electroweak loop corrections, dominated by the large 
negative mass shift from the top quark loop [31, 32, 33]. 

near-conformal resonance spectrum separated 
from light scalar 

moving up to 2-3 TeV with refined scale setting 

within reach of LHC Run 2 ? 

1 TeV 

light scalar at few hundred GeV? 

EW self-energy shift 
observed Higgs-like? 

dMjj ~ -12K^r]m] ~ 




A 

ai 
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Figure 5: Schematic view of the emerging resonance spectrum. The parameters K and n are defined in [31], 


10 



















Toward the minimal realization of a light composite Higgs 


Chik Him Wong and Julius Kuti 


6. Mixed action on the gradient flow and the epsilon regime 


The alternative to safe extrapolation from a regime of competing scalar and pion masses to 
the massless fermion limit requires crossover to the epsilon regime of ^PT at low enough scale 
A where Goldstone dynamics begins to decouple from the scalar state. This is difficult to do and 
requires significant resources. To control taste braking we cannot go to lattice spacings larger 
than the one set by j3 = 3.20. The value of F ~ 0.025 at this lattice spacing requires large 48^ x 96 
lattice volumes to control the F • L > 1 condition which is necessary for convergent expansion in all 
regimes of ;t;PT, including the epsilon regime. Even for our largest V = 48^ x 96 and V = 40^ x 80 
lattice volumes control with F-L ~ 1 — 1.2 is just barely sufficient. For the lowest fermion mass 
m = 0.002 the scaling variable mEV ~ 100 is very large characterizing the p-regime of ^PT we 
used earlier in the analysis. Reaching the epsilon regime would require an order of magnitude 
decrease in the scaling variable mZV which presents a considerable algorithmic challenge and 
requires substantial resources. Decreasing the fermion mass an order of magnitude to m = 0.0002 
would increase the cost, scaling with ~ 1/m and calling for algorithmic improvements. 


taste splitting of Q=-1 Dirac spectrum 
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Figure 6: The upper right panel demonstrates the correct index theorem and the lower right panel displays the lowest 
eigenvalues in the epsilon regime. 

We developed a promising new strategy to overcome the problem by performing ;t;PT analysis 
in the crossover to the epsilon regime with partial quenching and a related mixed action. We take 
the p-regime gauge configurations of the lowest fermion masses on the largest lattice volumes and 
analyze the fermion condensate and the Dirac spectrum after gradient flow times t = 2 or 3 with 
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the valence fermion action where the original gauge link variables are replaced with the ones at 
flow time t. This strategy can be viewed as a mixed action based analysis. The first encouraging 
results are shown in Figure 6. The upper left panel shows the infrared part of the Dirac spectrum 
on the original gauge configurations with strong taste breaking evidenced by the absence of de¬ 
generate quartets. After gradient flow time t = 3 degenerate eigenvalue quartets emerge with the 
correct count of the topology dependent zero modes from the index theorem showing restored taste 
symmetry in the valence action. The lower left panel illustrates the degeneracy of the Goldstone 
pion with non-Goldstone pions (scPion in the plot). The lower right panel shows that the epsilon 
regime is reached with the scaling variable AZflowV ~ 10 where the fermion mass is replaced by 
the scale of the Dirac spectrum (m —)■ A) and Tflow. not RG invariant, is reduced by almost a factor 
of ten. 


7. The scale dependent renormalized coupling and beta function 


An important and independent consistency condition of the model would be provided by 
matching the scale dependent renormalized coupling of the perturbative regime to the scale depen¬ 
dent coupling of the non-perturbative phase associated with ;fSB. We proposed a gauge coupling 
earlier g(/i = 1/L), running with the scale set by the finite volume [34] and defined on the gradient 
flow of the gauge field [35]. Since the gradient flow at flow time t probes the gauge field on the 
scale the running coupling can be defined as a function of L in finite volume V = while 
holding c = (8t)'/^/L fixed: ac(L) = 47r(t^E(t))/[3(l -|- 5(c)]. This volume dependent coupling 
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Figure 7: Two different schemes for scale dependent couplings are illustrated with explanation in the text. 
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is particularly suitable to study the perturbative regime and departures from it [34] . The measured 
renormalized couplings are very accurate and the scheme defines a one-parameter family which 
can be adjusted for different goals [36, 37]. For illustration, preliminary post-conference results of 
the step beta function for the sextet model with two flavors are shown in the lower right panel of 
Figure 7. 

In the non-perturbative phase with ;^SB we are interested in a scale-dependent and volume 
independent renormalized coupling. At fixed lattice size, bare coupling, and fermion mass m we 
determine the appropriate flow time t(g^,m) to match any targeted flow-dependent renormalized 
coupling g^ calculated from (t^E(t)). Assuming that the footprint of the operator on the gradient 
flow is sufficiently small compared to the Compton wavelength of the pion for p-regime analysis, 
the dependence of t(g^,m) on m can be replaced by t(g^,M|) in ;^PT of pion dynamics with linear 
dependence of t(g^,M|) on in leading order [38]. Any residual finite volume dependence can 
be corrected in ;tPT. 

This strategy is illustrated by the step by step procedure in Figure 7. The upper left panel shows 
the determination of the flow time t(g^,m) of the targeted coupling g^ and the upper right panel is 
in surprisingly good agreement with the linear behavior in M^. At two different lattice spacings 
in the m = 0 chiral limit, the lower left panel shows the scale dependent renormalized coupling 
g^(to) as a function of scale variation with to. A scale dependent and volume independent step beta 
function can be determined from this procedure. A more comprehensive analysis of the data is paid 
of our ongoing investigations including the extrapolation of the step function to vanishing lattice 
spacing and matching the two different scale dependent couplings of Figure 7. 
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